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The recent success of the NIST group [T\ f2] 
in generating abelian gauge field in cold atoms 
has created opportunities to simulate electronic 
transports in solids using atomic gases. Very re- 
cently, the NIST group has also announced in 
a DARPA Meeting the creation of non-abelian 
gauge fields in a pseudo spin-1/2 Bose gas. While 
there have been considerable theoretical activi- 
ties in synthetic gauge fields, non-abelian fields 
have not been generated until now [3-16J. Here, 
we show that in a non-abelian gauge field, a 
spinor condensate will develop a spontaneous 
stripe structure in each spin component, reflect- 
ing a ground state made up of two non-orthogonal 
dressed states with different momenta. Depend- 
ing on interactions, this ground state can re- 
duce back to a single dressed state. These mo- 
mentum carrying stripes are the macroscopic 
bosonic counterpart of the spin-orbit phenomena 
in fermions that are being actively studied in elec- 
tron physics today. 

The beauty of the NIST scheme is its simphcity. It 
also illustrates the important fact that both abelian and 
non-abelian fields can be related through a single fam- 
ily of hamiltonian, and suggests in our opinion a "gen- 
eralized adiabatic" scheme for generating gauge fields 
with increasing complexity. The scheme is the follow- 
ing. Consider the hamiltonian h = p^/2M + W{r) that 
operates on an atom with internal degrees of freedom, 
such as alkali atoms with hyperfine spin F. W is a spa- 
tially varying potential in spin space with typical wave- 
vector q. The energy scale for spatial variation of W 
is then Cq = h q^ /2M. If W has a group of L states 
(L < 2F + 1) at the bottom of its spectrum lying within 
an energy range /S.E <C e^ and is well separated from all 
other higher energy spin states by e^, then the low energy 
phenomena of the system can be described within this re- 
duced manifold of L states. By going into a frame in this 
manifold that transforms away the spatial variations of 
the spin states, a gauge field emerges. The gauge field 
is abelian if L = 1, and non-abelian if L > 2. Thus, by 
successively moving the high energy states across e^ into 
the low energy group, one can increase the dimensionality 
of the low energy manifold and create non-abelian gauge 
fields with increasingly rich structure. It should be noted 
that this is very different from the A-scheme in most the- 
oretical proposals, which makes use of a set of dark states 
sitting distinctly above a short living ground state of the 
system [5|. In contrast, the generalized adiabatic scheme 



makes use of the lowest energy states. It therefore elim- 
inates the collisional loss and hence intrinsic heating in 
the A-scheme. 

Before proceeding, it is useful to note the unique fea- 
tures of non-abelian gauge fields. In the abelian case, 
a constant vector potential has no physical effects as it 
can be gauged away completely. This is not true for the 
non-abelian case because of its non-commutativity As a 
result, a constant vector potential does matter. More- 
over, non-abelian gauge fields inevitably lead to spin- 
orbit coupling, so any potential (such as a confining trap) 
that alters particle trajectory also causes spin rotation. 
This immediately implies significant differences between 
bosons and fermions. For fermions, Pauli principle forces 
them into different orbital states. The spin structure of 
the system is then a result of the contributions of all dif- 
ferent occupied states. In contrast, bosons will search 
for or even construct (through interaction effects) an op- 
timum spin states, and magnify it to the macroscopic 
level. The current experiments at NIST have provided 
us opportunities to study macroscopic spin-orbit effects. 

(A) The NIST setup and the effective hamil- 
tonian: The NIST setup consists of two counter propa- 
gating lasers with frequency difference uo and momentum 
difference q^ directed along x toward a spin F = 1 Bose 
condensate of ^^Rb atoms. There is also a magnetic field 
directly along y with a field gradient. (See FigurefTl). The 
lasers induces a Raman transition in the atom, transfer- 
ring linear momentum q-k to the Bose gas while increasing 
the spin angular momentum by ?i at the same time. The 
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single particle hamiltonian is h{t) = ^ -\-W{t)^ where 
where F is the spin-1 operator. ?i\ is the quadratic Zee- 



man energy. TiVly 



hQn 



Gy is the Zeeman energy 



produce by the magnetic field along y: The Qq term 
is due to the constant magnetic field and the Gy term 
comes from the field gradient. Qr is the Rabi frequency 
in the Raman process. In the frame rotating in spin space 
along y with frequency cj, the hamiltonian becomes static 
H = h{t = 0)^ TiuF^, and is given hy H = p^ /2M + W, 
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and Qy = Qq — ^ -^ Gy. 



Eq.([2| shows that W has a very simple level structure 
in the frame in spin space rotating along y with angle 




FIG. 1. Schematics of the experimental setup at NIST. 
The Raman process consists of two lasers with wave vectors 
kop + qx and kop, frequencies uOop + cj and Uop impinging on 
the atomic cloud. Atoms excited by the laser will have their 
momenta increased by qx while their spin projection along y 
is changed by 1, as shown in the energy diagram at top. 



qx. For simplicity, let us take F = 1 and G = 0. The 
following cases are of particular interests: 
(i) Abelian case: This occurs when Qji ^ A, Cq/h^ with uj 
tuned close to Qq^ hence Qy ~ 0. The ground state in the 
rotating frame is the m = -\-l state along z, isolated from 
other two states (m = 0, — 1 along z) by ^ hftR > Cq. 
(ii) Non-abelian case: This occurs when A :^ Qr, Cq^ with 
(J tuned closed to cj = l^o — A. In this case, the states 
m = 1 and m = along y lie at the bottom of the 
spectrum, separated from the third state tti = — 1 by 
2hX > Cq. We shall from now on focus on this case. 

Let ipl^ and ^J^ be the operators that create a bo- 
son with spin projection m along y in the laboratory 
frame and in the rotating frame in spin space; and 
'^m = (e*^^^^(^) = e*^^^(^^. Focusing on the lowest 

two states 6m, i 



1,0, the hamiltonian is 
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where hm = 4>ln4>m^ ^ = Y^m^^^ ^ = \Muj'^y'^ is an 
harmonic trap, fi is the chemical potential, Qran are inter- 
actions between bosons in spin states m and n, qiq = ^oi, 
and 
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When G = 0, the Schrodinger equation 



(4) 
(5) 



has the following property. If x is solution of Eq.Q, then 



XL{x)=e^^e-^''^{n)mnX*ni^), ^=1^ 



1 




(6) 



is also a solution, where 7 is an arbitrary phase. 

(B) Single particle ground state: In zero field gra- 
dient G = 0, the momentum eigenstates is of the form 

Xm {x) = e'^^Xm, X = iv)^ and Eq.([5| becomes 
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where we have defined 

Q = q/2, k=p 
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and have expressed hftR in terms of a wave- vector £ and 
angle for later use. 
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The eigenvalues come in two branches, with energies 
^miP) = M (^ + i-)VMF+^)- (See Figure 



El). The ground states are the minima of -E'o(p) at 



p± = ±ko - q/2, ko = +v/Q2 - ^VQ^ = (9/2)cos^, 
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with energy 
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The energy of the upper branch at these momenta is 
Ei{p±) = eq- h'^£^/{2MQ'^), which is higher by Cq. It is 
worth noting that the value of the ground state energy is 
not of the order —TiVLr^ but a higher energy —(TiVtRY/eq. 
The wavefunctions at these degenerate ground states are 
X^^^\x) = e^^^-X^^^\ 
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Note that the states x^^^H^) ^^^ connected by Eq.(p| 
with 7 = 7r/2. They are orthogonal due to their differ- 
ent momenta. The spin states, however, have non-zero 
overlap, since 



(p+|p_)=x'^+)tx(P-)=sin0. 



(13) 



(C) Pseudo spin-1/2 spinor condensate: Con- 
densing in the dressed states |p^^^), the field operator, 
which has the expansion (j)rn{x) = Xl^Xm (^)^p, turns 
into a spinor field of the form 



<^^{x) = A+xl^+\x) + A_xl^-\x). 



(14) 
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FIG. 2. The energy levels Eo(p) and Ei{p) as a function of 
A; = p + q/2. The lower branch ^o(p) has two degenerate 
minima dX k — =b/co, where Uq — {q/2)cos0. The energy 
difference between the lower and upper branch at di/co is Cq — 
n^q^/2m. 



Because of the non-zero overlap Eq.(13), the density 
^m(^) = |^m(^)P of each spin component will develop a 
stripe structure. This can be seen by noting that the to- 
tal density n{x) = ni{x)-\-no{x) and the "magnetization" 
m{x) = ni{x) — no{x) are given by 



n{x) = \A+\ 



sini9(A;A_ 






m{x) = — cos^d^^ 



\A.\'). 



+ C.C.) (15) 
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Note also that m{x) is independent of 0. Eq.(15) shows 
that the contrast of the oscillation is set by the over- 
lap sin^, whereas the wavelength of the stripe is ir/ko = 
27r/(g'cos^). Thus, both contrast and wavelength increase 
with for < 7r/2. 

The amplitudes A± are determined by minimizing the 
Gross-Pitaevskii functional of Eq.([3|, which is obtained 
by replacing ^m(^) with the c- number $^(x), and nm{x) 
with nm{x) = \^m{x)\\ Defining |Ap = \A^\^ + |A_p, 
and a± = A±/\A\, the GP functional then reads, 



/C = (Eo-m)|A| 



|A|'*G(a+,a_ 



(17) 



where |A|^G(a+,a_) = J gmn^m{x)nn{x). Note that 
while the contributions of A+ and A- are separated out 
in the kinetic energy term due to the different momenta 
of X*^^^H^) ^^^ x''^~\x)^ they are coupled through in- 
teraction due to the overlap of their spin functions. For 
example, J n1{x) = J[n^{x)-\-m^{x)-\-2n{x)m{x)]/4:^ and 
the mixing of A^ and A- appears in j ii?{x). To mini- 
mize /C, we first minimize G{aj^^a-) with the constraint 
|a+p + |a_p = 1 to obtain the optimal value {a^j^^aZ_)] 
and then we have 




FIG. 3. The phase diagram of pseudo spin 1/2 Bose gas: Re- 
gion I is a superposition of two dressed state with momentum 
p+ and P-, II and III are the single dressed states p+ and 
P- respectively, a, /3, ttc, and /3c are defined in text. 



Since the minimization is straightforward, we shall 
only present the results, which are shown in Figure 3. 
The phase diagram depends on the parameters 

<^ = ^io/^, /^ = (^ii-^oo)/^, ^ = (^ii+^oo)/2. (19) 
and two numbers ac and /3gjderived from the laser pa- 
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cos6>(2 - tan^6>). For ^11,^00,^10 > 0, (as 
there are three possibilities: (I) Two dressed 



states, with both A± ^ 0; single dressed state with (II) 
Xp+(x), (A_ = 0), or (III) Xp-{x), (A+ = 0). 

Phase (I) occurs within the triangle shown in Fig. 3, 
bounded by the lines xyc^yxc = XcVc- The region exists 
only when ac > 0, which means sinO < y^2/3. Other- 
wise, interaction effect will drive the condensate into a 
single dressed state. In phase (I), the amplitudes are 
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and Go = G{al,al) = - 2(2-2a-(i+c.)tan2 ^) +(1 + ^)(1- 
^ sin^ 0) . The relative phase between A+ and A- , how- 
ever, remains undermined within the GP approach. This 
phase can be fixed by perturbations such as field gradient 
the breaks the symmetry Eq.(p|, or by quantum fluctu- 
ation effects that go beyond GP approach. As discussed 
before, the density of each of the spin component rii and 
no of this phase has a stripe structure. The case /3 = 
(^11 = ^00) is special. In that case, we have \A_^\ = \A_\ 
for a < ac. For a > ac, the two dressed states x^^^^ ^^^ 



are degenerate. 
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Ap- 
In the presence of a harmonic trap V{t) = ^Muj^t 
|A|^ = (/i - Eo)/Go, Go = G{a''_^,a''_). (18) with harmonic length d = y^h/{Muj) > 27T/q, the wave- 



length of the stripe, we can apply Thomas-Fermi ap- 
proximation, and the condensate wavefunction is given 
by Eq.(14), (18) and ( [2Q| ) with chemical potential /i in 
Eq.(18) replaced by /i(r) = /i — V{r)^ i.e. for Phase (I), 



text. Let us return to the single particle hamiltonian 
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The density profile ni(r) for the m = 1 spin component 
along y is shown in Figure 4 for the parameters = ^tt 
AT = 2.5 X 10^ etc. 

Apart from the stripe structure, the presence of these 
phases can be detected by measuring the displacement 
of the atomic cloud in expansion after the trap is turned 
off. For single dressed state with momentum p+ and p- , 
the atom cloud will displace in x-direction by an amount 
determined by p+ or p_. For the condensate with two 
dressed states, the system after expansion will separate 
into two atom clouds moving with different momenta. 
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FIG. 4. The upper figure is the column density ni{x,y) = 
J dzni{x,y, z). The lower frame is ni(x,0). The period of 
oscillation is n/ko. The contrast of oscillation at the center is 
70%. Our calculation is performed for ^^Rb with iV = 2.5 x 



10^ atoms, q = 1.56 x Wm~ 



hQf 



h X T.lKHz, 



cloud size Rtf = 20/im 1 . The values a and /3 used are given 
by a = |ac and /3 = |/3c. The length displayed is in units of 
the laser wavelength 804. 3nm in ref.^lj . 

Our discussions here focus on current experiments. It 
is, however, useful to put them in a more general con- 
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= S, 



^nP^/2M + Wmn(r), where W(r) is a matrix 
in a spin-F space, m = — F, —F + 1, • • • , F. Let |n)r be 
the eigenstate of W with energy Sn^ which we label as 
Sf < £f-i < ••• < ^-F in ascending order. The unitary 
matrix U{r) that diagonalizes W in the original spin basis 

is Umn{T^) = {m\n)r, and [U"^ {t)W {t)U {T)]mn = SmSmn- 

In the frame where W is diagonal, the hamiltonian is 

[U{T)HW{T)]mn = 2!^ [(P + A)^]^^ + ^n^mn, whcrC the 

Amn = (^^|rP|^)r IS the gaugc field emerging in this 
frame [17 . The gauge field is non-abelian if [A*, A^] ^ 0. 
For the NIST experiment, one restricts to the lowest two 
states of spin-1 system. When the field gradient G 7^ 0, 
A has both x and y component and is non-abelian. The 
G = case that we considered is special, since A has 
only one component. However, this uniform gauge field 
can not gauged away without changing the hamiltonian 
due to the presence of the Rabi term Vtp>T\ in Eq.Q. 
The phenomena associated with the non-commutativity 
A caused by field gradients are sufficiently rich that they 
deserve separate discussions. 
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